Clustered chimera states in delay coupled oscillator systems 
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We investigate chimera states in a ring of identical phase oscillators coupled in a time-delayed 
and spatially non-local fashion. We find novel clustered chimera states that have spatially dis- 
tributed phase coherence separated by incoherence with adjacent coherent regions in anti-phase. 
The existence of such time-delay induced phase clustering is further supported through solutions 
of a generalized functional self-consistency equation of the mean field. Our results highlight an 
additional mechanism for cluster formation that may find wider practical applications. 
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The study of time delay induced modifications in the 
collective behaviour of systems of coupled nonlinear 
oscillators is a topic of much current interest both for its 
fundamental significance from a dynamical systems point 
of view and for its practical relevance to modeling of 
various physical, biological and chemical systems. In real 
life situations time delay is usually associated with finite 
propagation velocities of information signals, latency 
times of neuronal excitations, finite reaction times of 
chemicals etc. and in collective oscillator studies it is usu- 
ally modeled through a time delayed coupling function. 
Many such past model studies on globally and locally 
coupled oscillator systems have uncovered interesting 
and sometimes novel time delay induced modifications 
of the equilibrium, stability and bifurcation properties 
of their collective states [l[ . While global and local 
(nearest neighbour) coupling models have traditionally 
received much attention there is now a growing interest 
in the collective dynamics of models with non-local 
couplings [3, d, 0, 0, @, 0]. Non-local coupling can 
be relevant to a variety of applications such as in the 
modeling of Josephson junction arrays Q, chemical 
oscillators 0, S, , neural networks for producing snail 
shell patterns and ocular dominance stripes [13, EH, E3] 
etc. One of the striking features of non-locally coupled 
oscillator systems is that they can support an unusual 
collective state in which the oscillators separate into two 
groups - one that is synchronized and phase locked and 
the other desynchronized and incoherent Such a 

state of co-existence of coherence and incoherence does 
not occur in cither globally or locally coupled systems 
and has been named as a chimera state by Abrams & 
Strogatz [3]. The nature and properties of this exotic 
collective state as well as its potential applications are 
still not fully explored or understood and therefore 
continue to offer exciting future possibilities. It is not 
known for example whether such chimera states can 
exist in the presence of time delay in the system and 
if so then what their characteristics are. This is the 
principal question we examine in this work through 



numerical simulations and mathematical analysis of 
a model system consisting of a ring of densely and 
uniformly distributed identical phase oscillators that 
are coupled in a time-delayed and spatially non-local 
fashion. We find that chimera states do indeed exist 
but acquire an additional spatial modulation such that 
the single spatially connected phase coherent region of 
the usual chimera state is now replaced by a number 
of spatially disconnected regions of coherence with 
intervening regions of incoherence. Furthermore the 
adjacent coherent regions of this clustered chimera state 
are found to be in anti-phase relation with respect to 
each other. To understand the origin and the nature of 
this pattern we have extended the mean field approach 
used by Kuramoto Q and applied it to our system 
which has a distance-dependent time delay factor in the 
coupling and have derived a functional self-consistency 
equation. A numerical solution of this self-consistency 
equation yields a space-dependent order parameter and 
a space-dependent mean phase function that confirm the 
existence and explain the nature of the spatial pattern 
of the oscillator phases. 

We consider the following model equation representing 
the continuum limit of a chain of identical phase oscilla- 
tors arranged on a circular ring C, 

(x, t) = u — G(x — x) 
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x sin[4>(x, t) — <j)(x , t — t x x i) + a] dx' (1) 



where 2L is the system length, u> is the natural frequency 
of the oscillator and a closed chain configuration is en- 
sured by imposing periodic boundary conditions. The 
kernel G{x — x'), appropriately normalized to unity over 
the system length, is taken as, 

Gix-x^^J-^e^' (2) 

which provides a non-local coupling among the oscilla- 
tors over a finite spatial range of the order of k^ 1 which 
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is taken to be less than the system size. The coupling is 
time delayed through the argument of the sinusoidal in- 
teraction function, namely, the phase difference between 
two oscillators located at x and x' is calculated by taking 
into account the temporal delay for the interaction signal 
to travel the intervening geodesic (i.e. shortest) distance 
determined as d x ^ x i = min{\x — x'\, 2L — \x — x'\}. The 
time delay term is therefore taken to be of the form, 
t x ,x> = d x ,x> /v where v is the signal propagation speed. 
In the absence of time delay the above equation reduces 
to the one investigated in [5j, [6| . The constant phase shift 
term a in the undclayed model breaks the odd symme- 
try of the sinusoidal coupling function and as discussed 
in (l3l . 14 1 it is needed as a tuning parameter for ob- 
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taining chimera solutions in the undelayed case. In the 
presence of time delay however we find that a no longer 
plays such a critical role since the time delay factor also 
fulfills a similar function. 

We now describe direct numerical simulation results 
obtained by solving Eq.JI]) using a large number of dis- 
crete oscillators (typically N=256). The set of system 
parameters chosen for the simulations illustrated here 
were, 2L = 1.0, a = 0.9, fc" 1 = 0.25, lo = 1.1 and 
v = 0.09765625 corresponding to a maximum delay time 
(Tmax) in the system of 5.12. As discussed in past stud- 
ies [D, @], the choice of appropriate initial conditions is 
very important for numerically accessing a chimera state. 
Kuramoto used a random distribution with a Gaussian 
envelope for the initial distribution of the phases to ob- 
tain a chimera solution. For our time-delayed system 
we find that choosing the initial phases of the oscillators 
from a uniform random distribution between and 2n 
and then arranging them in a mirror symmetric distri- 
bution in space provides a rapid access to a clustered 
chimera state. Our simulations have been done with the 
XPPAUT [l5| package using a Runge-Kutta solver (with 
a small integration time step of St = 0.01) till a time 
stationary solution is obtained and tested for indepen- 
dence from discreteness effects by repeating the runs for 
N = 128, 256 and 512. In panels (a) and (b) of FigQ] 
we show a space time plot of our simulation for the pa- 
rameters mentioned above in the early stages of evolu- 
tion (starting from random initial phases) and in the fi- 
nal stages of the formation of a clustered chimera state 
respectively. Panel (c) shows a snapshot of the spatial 
distribution of the phases in the final stationary state. 
We see four coherent regions interspersed by incoherence 
and also note that the adjacent coherent regions are in 
anti-phase. Panel (d) is a blowup of the region between 
x = —0.5 to x = —0.25 giving an enlarged view of an 
incoherent region and portions of the adjacent coherent 
regions. These solutions are also found to be quite robust 
and show no signs of instability over arbitrarily large in- 
tegration times. We have tested the integrity of the solu- 
tions for times well over l00T max . A detailed parametric 
study of the stability regions is presently under progress 




FIG. 1: (color online) (a) The space-time plot of the oscillator 
phases cj> for the parameters 2L = 1.0, k = 4.0, 1/v — 10.24, 
uj = 1.1 and a — 0.9 in the early stages of evolution from 
a random set of initial phases. Panel (b) shows a later time 
evolution and panel (c) gives a snapshot of the final stationary 
state. Panel (d) is a blowup of the region between x = —0.5 
to x = —0.25 giving an enlarged view of an incoherent region 
and portions of the adjacent coherent regions. 



and will be reported elsewhere. 

To gain a better understanding of the nature of this 
pattern and of the dynamics of its formation we have 
carried out a mathematical analysis based on the gener- 
alized mean field concept as developed by Kuramoto for 
the non-delayed case. For this we first rewrite Eq.(TT]) in 
terms of a relative phase 0(x, t) = <p(x, t) — Qt (where Q 
represents a rotating frame in which the dynamics simpli- 
fies as much as possible such that with the phase-locked 
portions rotate with this constant drift frequency) as, 

^e(x,t)=u-Q- J G(x-x') 

x sin[0(x, t) — 9(x',t — t x ,x') + ct + ilr x ,x'] dx' (3) 

The key idea behind Kuramoto's analysis of chimera 
states was the introduction of a mean field like quan- 
tity, namely, a complex order parameter Re 1 ® , defined in 
a manner analogous to what is done for globally coupled 
systems. For our case we write, 

R(x,t)e ie ^= J G(x - x'yW't-^-^'Ux' 

(4) 

The above order parameter differs from the usual defi- 
nition for global coupling systems in several ways - the 
spatial average of e %B is weighted by the coupling kernel 
G(x — x 1 ), the phase 9 is evaluated in a time delayed 
fashion and the factor e _ *° Tx x ' adds a complex phase to 
the kernel G(x — x'). The latter two features provide a 
further generalization of Kuramoto's analysis carried out 
for a non-delayed system @, @, [l3, 14 1. 
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In terms of R and 0, Eq.© can be rewritten as 
d 

— 6(x, t) = A — R(x, t) sin[0(x, t) - 0(x, t) + a] 
at 



(5) 



where A = u> — 0. Eq.© is in the form of a single 
phase oscillator equation being driven by a force term 
which in this case is the mean field force. To obtain a 
stationary pattern (in a statistical sense) we require R 
and to depend only on space and be independent of 
time. Under such a circumstance the oscillator popu- 
lation can be divided into two classes: those which are 
located such that R(x) > A can approach a fixed point 
solution (dO(x,t)/dt = 0) and the other oscillators that 
have R(x) < |A| would not be able to attain such an 
equilibrium solution. The oscillators approaching a fixed 
point in the rotating frame would have phase coherent 
oscillations at frequency f2 in the original frame whereas 
the other set of oscillators would drift around the phase 
circle and form the incoherent part. Following the pre- 
scription provided by Kuramoto [B], Q for the undelayed 
case, we substitute the solutions of Eq.© for the two 
classes of oscillators into the integrand on the right hand 
side (R.H.S.) of Eq.© and obtain the following func- 
tional self-consistency condition, 



R{x)e i@ ^ = e if} j G(x - x 'yl^')-^, x> ] 



A-^WW) 
R(x') 



(6) 



where f3 = ir/2 — a. We need to solve for three unknowns 
- the functions R(x), 0(x) and the quantity A. Con- 
dition © provides only two equations when we separate 
its real and imaginary parts. A third condition can be 
obtained by exploiting the fact that the equation is in- 
variant under any rigid rotation <d(x) — > 0(x) + Go- We 
can therefore specify the value of 0(x) at any arbitrary 
chosen point, e.g. Q(L) = 0. We have solved Eq.© 
numerically by following a three step iterative procedure 
consisting of the following steps. We choose arbitrary but 
well behaved initial guess functions for R(x) and 0(x) 
and use the condition 0(£) = in one of the equations 
of © to obtain a value for A. The initial profiles and the 
A value so obtained are used to evaluate the R.H.S. of 
© to generate new profiles for R and 0. These are next 
used to generate a new value of A and the procedure is 
repeated until a convergence in the value of A and the 
functions R and are obtained. 

A MATHEMATICA program incorporating this algo- 
rithm was developed and benchmarked against the re- 
sults for the no-delay case. Fig[2]shows the rapid and ex- 
cellent convergence in A to a unique value of A = 0.189 
for the solution of Eq.© with system parameters cho- 
sen identical to the ones that were used to obtain a 
clustered chimera state by a direct solution of Eq.©. 
The converged spatial profiles of the order parameter (i? 
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FIG. 2: Variation of A with the iteration number showing 
a rapid convergence in the numerical solution of the self- 
consistency Eq.©. The system parameters are identical to 
those used in the direct solution of Eq.([lJ. 



and 0) are shown in Fig [3] and the converged value of 
A is marked in the upper panel by the horizontal line. 
The amplitude of the order parameter (R) shows a peri- 
odic spatial modulation - peaking at four symmetrically 
placed spatial locations. The corresponding phases of 
the order parameter are seen to be in anti-phase for ad- 
jacent peaks in R. In between the peaks R is seen to 
dip to very small values at certain locations such that 
R(x) < |A| which should correspond to the incoherent 
drifting parts of the chimera. To better appreciate the 
agreement between the direct solutions of Eq. © and the 
mean field solutions of Eq. ([6]) we have plotted the results 
together in FigJU As is clearly seen the measured order 
parameter (R and 0) and A from the direct simulations 
of Eq.© match well with the results of solving Eq.©. 
The spatial profile of the phases (<f>) of the oscillators as 
obtained from the direct simulation of Eq.© is shown 
in the top panel of Fig[U We see four coherent regions 
interspersed by incoherence as expected from the results 
of solving Eq.©. 

We note from Fig.© and Fig.© that both R and 
are mirror symmetric (i.e. R(x) = R(— x),0(x) = 
0(— x)), a property that the original phase Eq.© also 
possesses. Eq.© is also invariant under the transforma- 
tion (</>(x, t) — > —(f)(x, t), uj — > — u>, a — » —a) and can have 
solutions with such a symmetry as well, namely, travel- 
ing wave solutions given by 0(x,t) = Ctt + irqx/L. In 
our numerical simulations we find that by changing the 
initial conditions, but keeping the same system parame- 
ters, we can also get traveling wave solutions. There also 
seems to be a clear correspondence between the number 
of clusters of the observed chimera state and the wave 
number q of the co-existent traveling wave solution. For 
the 4-cluster chimera of Fig|3] the co-existent traveling 
wave has q = 2 and similar results have been obtained 
for 6-cluster (q = 3) and 8-cluster (q = 4) chimera solu- 
tions. 
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FIG. 3: Spatial profiles of the amplitude R and the phase O of 
the order parameter obtained by solving the self-consistency 
Eq.© by an iterative scheme. The horizontal line in the 
upper panel, marks the converged value A = 0.189. 



the resultant pattern of phase distribution seen in the 
clustered chimera state. Thus time delay offers an ad- 
ditional mechanism for cluster formations in dynamical 
systems and model systems incorporating time delay may 
provide a useful paradigm for studying this phenomenon. 
Our results can be usefully extended to higher dimen- 
sions e.g. to examine the influence of time delay on spi- 
ral wave based chimeras in two dimensions Q and may 
also help provide insights into experimental observations 
of clustered states that are generic to many chemical and 
biological systems [l^. 
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FIG. 4: (color online) (a) The phase pattern for a clustered 
chimera state as obtained by direct simulation of Eq.flJ. The 
measured spatial profiles of the order parameter (R and 0) 
from these simulations are shown in panels (b) and (c) as 
dashed curves and compared with the solutions from the self- 
consistency Eq.© shown as solid curves, (d) u> — cj> for the 
oscillators from a direct simulation of Eq.((l|. The horizontal 
lines in (b) and (d) mark the converged value of A = 0.189. 



To conclude, we have demonstrated for the first time 
the existence of chimera type solutions in a time-delayed 
system of non-locally coupled identical phase oscillators. 
Time delay is found to lead to novel clustered states with 
a number of spatially disconnected regions of coherence 
with intervening regions of incoherence. The adjacent co- 
herent regions of this clustered chimera state are found 
to be in anti-phase relation with respect to each other. 
Our numerical simulations are further validated and ex- 
plained through solutions of a generalized functional self- 
consistency equation of the mean field. The mean field 
parameters (the amplitude and phase of the complex or- 
der parameter) clearly reflect the modulated nature of 
the effective driving force on each oscillator and lead to 
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